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Abstract 

We explain the meaning of dynamical manipulation, and we illus- 
trate its mechanism by using a system composed of a charged particle in 
a Penning trap. It is shown that by means of appropriate electric shocks 
(delta-like pulses) applied to the trap walls one can induce the squeezing 
transformation. The geometric phases associated to some cyclic evolu- 
tions, induced either by the standard fields of the Penning trap or by the 
superposition of these plus a rotating magnetic field, are analysed. 

1 Introduction 

The quantum control, or wavepacket engeenering, is one of suggestive subjects 
in Quantum Mechanics (QM). The very name reflects the dreams of almost 
any physicist: to be able to induce a physical system to make anything we wish 
Since long time ago we have been involved with that subject using 
a specific name: dynamical manipulation problem. My goals in this work are 
diverse: 

• To explain what we mean precisely by dynamical manipulation 

• To illustrate how it works in a realistic physical arrangement, namely the 
Penning trap 

• To show the relation between dynamical manipulation and some 'hot' 
subjects in QM as squeezing and geometric phases 

x Talk given at the workshop 'Symmetries in quantum mechanics and quantum optics', 
Burgos (Spain), September 1998. To be published in the Proceedings. 
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Having this in mind, this work has been organized as follows. In section 2 some 
generalities of the manipulation problem will be presented. Some specific dy- 
namical processes, the evolution loops, will be introduced, and we will explain 
that they can be used in order to generate an arbitrary unitary transformation. 
In particular, they will be useful to induce the free evolution going back in time. 
In section 3 those general techniques will be applied to a realistic arrangement, 
namely, the Penning trap. Thus, it will be shown that the evolution loops occur 
in the Penning trap, and they will be called Penning loops. The 'perturbed' 
Penning loops will be used as a starting point to induce the squeezing transfor- 
mation (in general the scale operation) and the Fourier-like transformation. In 
section 4 the general setting of the geometric (Aharonov-Anandan) phase will 
be shortly presented. Section 5 will stablish the connection between the geo- 
metric phase and the Penning loops (perturbed and unperturbed). The paper 
will end with some general conclusions. 



2 Manipulation problem 

Motivation: at first sight there is an asymmetry in nature so that the evolution 
forward in time is privileged. One might think that certain unitary operators 
cannot be dynamically achieved, e.g., the inverted free evolution towards the 
past (backwards in time) 

e iX ^, A>0. 

The dynamical manipulation problem tries precisely to answer the question: 
can any given unitary transformation U on a physical system be dynamically 
induced? In other words, can we find a Hamiltonian H(t) such that U arises 
from a solution to the Schrodinger equation with such a H(t) at some time 

t — T? 

U = U(t = r), —U(t) = -iH(t)U(t), U(t = 0)=I. 
at 

The essence of the idea was formulated by Lamb in terms of system states, 
namely, given any two state vectors a) and |6) one looks for a Hamiltonian H(t) 
which links them dynamically so that there is a solution to Schrodinger 

equation becoming \a) and \b) at two different times t a , tb [Q. Later on, this 
idea was pursued at the operator level (as posed above) by Lubkin, Miclnik, 
Waniewski, the present author and other colleagues p[-[f[3"). 

There are signs that theoretically any unitary operator can be dynamically 
induced if two assumptions are made || : 

• Any potential V(q, t) continuous in (q, t) and leading to a self-adjoint 
Hamiltonian 

2 

H(t) = ^- + V(q,t) 
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can in principle be created; in such a case it is said that 



(1) 



where T is the time ordering symbol, are dynamically achievable operators 
(DAO) 

• The limit of operators of the form (1) are DAO 
Some simple examples of DAO are: 



• exp(— it^r), t > 0: induced by the free particle Hamiltonian H = 

• exp[— iV(q)]: induced by the Hamiltonian associated to the 'kick' of po- 
tential V(q, t) = S(t — to)V(q); it can be seen as the limit when e — > of 
the operator sequence 

U e = U(t , t + e) = e -H4+in«)] = e -He4+v( q )] . 



2.1 Evolution loops 

The evolution loops (EL) are specific dynamical processes such that the evo- 
lution operator of the system becomes / (modulo phase) at some time r 

U{t) = e^I = I. 

They are a natural generalization of what happens for the harmonic oscilla- 
tor. However, the interest on those processes arose after the discovery of the 
following operator identity ||] : 



-i*- e -ix-\ =/, A>0. (2) 

Notice that the left hand side of (2) involves just DAO so that the EL are also 
DAO. A schematic representation of (2) is shown in figure (l.a), where the 
sides of the figure represent intervals of the free evolution of length A while the 
vertices represent kicks of potential q 2 /2 of intensity 1/A. It is easy to show 
(2) in the Heisenberg picture, by considering the LHS as an evolution operator 
and noticing that 

c A 2 e 2 l ie 2 e * 2 ! ' 1 ' J 



o iA4 iAp 
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Figure 1: a) schematic representation of the evolution loop (2); b) the precession of 
a generic evolution loop induced by some perturbation 

2.2 Perturbed evolution loops 

The EL are important by the suggestive idea that by applying perturbations, 
the complete Hamiltonian (the loop Hamiltonian plus the perturbation) will 
induce the precession of the distorted loop which can lead, in principle, to any 
unitary operator This process is schematically represented in figure (l.b), 
where the gray line represents the EL and the black line the precession induced 
by the perturbation. 

Some additional interesting applications of the EL can be found. An impor- 
tant one arises from equation (2) by noticing that the free evolution backwards 
in time is also a DAO 



Another suggestive application concerns measurements in QM. To test, e.g., 
the reduction axiom one has to make basically an entire sequence of almost 
'simultaneous' measurements (so that the subsequent evolution will not destroy 
the wavepacket resulting of the first measurement). For a system in an EL, 
once the first measurement has been made we can assure that the reduced 
wavepacket, independently of what it was, will be reconstructed at the finite 
loop time r. So, we get more freedom to perform the second measurement, and 
the EL is a device avoiding that the reduced wavepacket will be demolished by 
the natural evolution [fL5|. 




A > 0. 



(3) 
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3 Manipulation and Penning trap 

Trying to find a realistic system in which the above techniques could be applied 
we arrived to the Penning trap. A charged particle inside an ideal hyperbolic 
Penning trap is under the action of a constant homogeneous magnetic field 
along z-direction plus an electrostatic field induced by a quadrupolar potential, 
characterized by the Hamiltonian || [T(| : 

# =^(p--A(r)) 2 + eV(r), 
2m v c / 

A(r) = -irxB , V(r) = V (r 2 - 3z 2 ). 

The region where the charge is trapped is limited by the equipotential surfaces 
of V(r). In the real trap, electrodes with that form are placed at the right 
positions (see our computer construction in figure 2). The two endcaps are at 
the same potential while the ring is at a different one. 




Figure 2: simulation of the Penning trap cavity. The missing section of the ring was 
removed to see better the inside of the cavity. 
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The characteristics of the motion can be easily seen by expanding explicitly 
the previous Hamiltonian: 

H = H z + H c + Hp, 
where H z , H c and H p are commuting Hamiltonians given by: 

H — 12 I to, .2 „2 

tf — " c r 
n c — y Zi 

^ = 2^ (f£ + Py) + f ^ (* 2 + 2/ 2 ) = mti + 
and the trap regime is guaranteed if 

o 4eVn e_B n o wj? — 2^n 

= " > 0, = — > 0, ^ = c > 0. 

m mc 4 

Along z direction we have harmonic oscillator motion of frequency loq. On the 
x — y plane the motion consists of rotations around z axis with frequency uj c /2 
superposed to a harmonic oscillator motion of frequency tu p along p. Notice 
that there are indeed just two independent frequencies. 



3.1 Penning loops 

We want not only trapping, but to induce some finer manipulations on the 
system, e.g., the evolution loops. They arise if the three motions H z , H c , H p 
are 'sinchronized', i.e., if the frequencies u)q, lu c , lo p are commensurable. There 
are different possibilities: 

• If oj p /ljq — I2/W £ Q, the z — p motions can be sinchronized. What is 
left at r = l\T, T = 2tt / loq are pure rigid rotations 

• If uj c /(jJq G Q, it is not guaranteed that lj p /ljq will be also rational. 
However, there are some values of u> c /loq G Q for which this happens 

w c /w = 3/2, uj p /uj = 1/A, t = 2T, 
lo c /u = 9/4, u p /u = 7/8, t = AT, 
ll> c /luo — 33/8, u) p /u>o — 31/16, r = 8T. 

The results above mean that the evolution loops are DAO of the charge inside 
the Penning trap. We have called them Penning loops (PL) H. 

3.2 Perturbed Penning loops 

From now on, let us restrict ourselves to the PL with period r = 2T, i.e., 
take lo c = 3a->o/2, uj p = u)q/A. Let us 'perturb' this PL by a sucession of two 
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instantaneous discharges applied to the walls of the trap, represented by the 
potential 

V"(r, t) = [Vj6(t - ti) + V^5(t - t 2 )\ (r 2 - 3z 2 ) h < t 2 < r. 
The total Hamiltonian is again a sum of 3 commuting terms 
H(t) = H l + eV'(r,t) = H z (t) + H l c + H„(t), 

H z (t) =H l z + ^ [F'S(t - h) + F"8{t - t 2 )] z 2 , 
H p (t) = ffp - j [F'<J(t - ti) + F"<5(i - t 2 )] P \ 

TO ' TO ' 

where the superindex I means to choose the specific values of loq, ui c and to p 
producing the EL with r = 2T. The evolution operator at r takes the factorized 
form 

U(r) = U z (t)U c (t)U p (t), 

where 

U z (t) = e -iH l jT-t2) e -iF''mz 2 /2 e -iH l Jt 2 ~t 1 ) e -iF'mz 2 /2 e -iH l z t 1 ^ 

U c (t)= e- iH 'c T = e i3vL * , 

U p (t) = e" iff P (r "* 2) e lF " mp2 / 4 e - lff p( t 2-*i) e iF'mp 2 /A e -iH l p ti _ 

For quadratic Hamiltonians, J7(r) is defined by the linear transformation 
^(r)(p)^) = u(r)(p), 

where u(t) is a 6 x 6 simplectic matrix called evolution matrix. In our case 
that matrix reduces to 2 x 2 unimodular matrices u x = u y = u x (t) = u y (r), 
u z = u z (t) acting on the pairs {x,p x ) T ', (y,p y ) T , (z,p z ) T : 

u x = -u ho {±f,T - t 2 )u k {-^-)u ho {±f,t 2 - t 1 ) Ufc (-^l) Ufco (SflL ) tl ), 

(4) 

u z = Mho(wo, t - t 2 )u k {F")u ho (uo,t 2 - h)Uk(F')Uho(U0, h). 

Notice that the evolution matrices U} Wl u k , and the corresponding operators 
inducing them are of the form: 

(r.r^a , ,-f sinqjf \ 2 „ 2 

<*»«* — \ (i)=e _ l(fcW¥)t 
-ratjsmwt coswt y w 

MF)=(_n F j)> tr fc (t ) = e- im ^. 

We are interested in the following transformations (the corresponding pa- 
rameters making them true are immediately reported in the corresponding 
table) : 
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• 3-dim Fourier-like transformation 

x \ 2 p x y -> Mp y z 
V* -» —&x p y -» p z 









TTIUJq A i 


mw A 2 


5.3131 


7.2533 


1.5165 


3.5238 


-39.0332 


5.3131 


7.2533 


-1.5165 


-0.6046 


6.6970 


0.9701 


11.5962 


1.5165 


0.6046 


-2.3891 


0.9701 


11.5962 


-1.5165 


-3.5238 


0.4099 



Table 1. The parameters producing a 3-dim Fourier- like transformation. 

• 1-dini Fourier-like transformation in z plus 2-dim scale transformation in 
x-y 

x — > A 2 .t y — > \ 2 y z — > Aip 2 

Px -> yT/Px ~ » JjPy Pz ~ > ^XT^ 



w ii 


^0*2 


F /loq — F /loq 


mwoAi 


A 2 


1.2094 


5.0738 


-2.1381 


-6.3874 


-10.2781 


1.9322 


4.3510 


-2.1381 


-1.0959 


-3.6353 


7.4926 


11.3569 


-2.1381 


-6.3874 


-0.0973 


8.2153 


10.6342 


-2.1381 


-1.0959 


-0.2751 



Table 2. The conditions to produce 1-dim Fourier-like and 2-dim scale trans- 
formations on z and x — y respectively. 

• 3-dim scale transformation 

x — > \ 2 x y — ► A 2 y 2 — ► Aiz 

Px-»^Px Pv-£p» Pz^xtp* 





W0<2 


F'/uq 


F"/wo 


Ai 


A 2 


1.2363 


8.4896 


-1.1589 


0.7524 


0.4712 


5.0901 


2.2064 


7.5194 


-0.7524 


1.1589 


2.1222 


5.0901 


4.0768 


11.3301 


0.7524 


-1.1589 


2.1222 


0.1965 


5.0469 


10.3600 


1.1589 


-0.7524 


0.4712 


0.1965 



Table 3. The parameters inducing the 3-dim scale transformation. 



To find the values producing, e.g., the 3-dimensional Fourier- like transfor- 
mation, we impose that the off diagonal elements of u x and u z of (4) become 
null. By solving these equations we will find the values of u>oh, i^oh, F'/cjq 
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and F" /luq inducing that transformation, and from these it is simple to find 
the corresponding values of the parameters Ai and A2 (see table 1). The same 
procedure can be used for the other two transformations (see tables 2 and 3) 



We are specially interested in the scale transformation due to its close con- 
nection with squeezing. From table 3, it is clear that we can produce different 
combinations of squeezing and amplification. For instance, if we take the values 
in the first row, at the end of the full process we will have produced a squeezing 
with scaling 0.47 along z direction and an expansion on x — y scaled by 5.09. 
By taking the values in the fourth row we will have produced the 3-dimensional 
squeezing with scalings 0.47 along z and 0.19 on x — y plane. It is interesting 
to notice that a particular scale transformation can also be gotten if 



F" F' woti 

— = = cot(— — ). 

uj uj 2 



oj t 2 = uJoti + 2n, 
The scaling parameters are in this case (see figure 3): 

/l+COs(^Lk)\ 2 



This means that along z direction an EL is again produced but on the x — y 
plane we have gotten the scale transformation. The squeezing arises when ui^ti 
takes values in the interval (n, 2n) with the corresponding kick intensities as 
given above. 



100 




Figure 3: the scale parameter A2 as a function of £ = Uoti. 
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4 Geometric phase 

The geometric phase was discovered by Berry, who realized that for cyclic adi- 
abatic evolution of the eigenstates of a slowly changing cyclic Hamiltonian, 
H(t) = H(Q), there is associated a geometric phase factor unnoticed by many 
people that had been working by years with the adiabatic approximation fl^| . 
Later on, Aharonov and Anandan realized that the key point of Berry's ap- 
proach was the cyclicity of the state rather than the adiabatic assumption, 
and thus they associated a geometric phase to any cyclic evolution regardless 
whether or not the Hamiltonian inducing the evolution changes slowly in time 
or even if it is time-independent [|18[ . The most economic formulation of the 
geometric phase runs as follows J_9| . 

Suppose that a system state \^p(t)) is cyclic, i.e. 

\i>(r)) = e^|V(0)>. 

It turns out that c/i is a sum of a dynamic plus a geometric contribution, the 
last one called geometric phase is given by 

p = ^+i J\m\j t m))dt. 

(3 is geometric in the sense that it is the holonomy of the horizontal lifting of 
the closed trajectory in the projective Hilbert space V, which arises due to the 
curvature of V . Hence, (3 measures global curvature effects of the projective 
Hilbert space. It is worth to notice that the calculation of the geometric phase, 
even the determination of the cyclic states of a given system, is not an easy 
task. This is one of the reasons why a lot of people have became involved in 
the subject H|, 0, 111111,11. 



5 Geometric phase and Penning loops 

As pointed above, for a system in a Penning loop any state \ip(t)) is cyclic with 
period equal to the loop period: 

|V(r)>=e*|^(0)>. 

Thus, it is quite natural that a geometric phase factor should be associated to 
any state. For evolution loops induced by time-independent Hamiltonians, the 
geometric phase is directly related to the expected value of the energy in the 
initial state: 

= ^ + T(ip(0)\H l \ip(0)). 

An alternative formula arises working in the basis of energy eigenstates of H 1 , 
where again the superindex / means that we are taking the parameters of H so 
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as to induce the corresponding Penning loop. Let be \E nin2n3 ) the eigenstate 
of H l associated to the eigenvalue E nin2n3 : 

H \E niTl2n3 ') E nin2Tl3 \E ni7l2n3 ) . 

By decomposing now \ip(0)} in that basis 



IV'(O)) — ^ ] C nirl2n3 \E niTl2rl3 ), 
nin 2 n 3 =0 



we finally get 



P = 4> + J2 E ^ 



li j 113 c ™i«2«3 • 



nin2«3=0 



This formula is similar to the one obtained for the evolution loop induced by 
the harmonic oscillator Hamiltonian |T3] (see also pi). 



5.1 Geometric phase and perturbed Penning loops 

Now, instead of 'perturbing' the PL by a term affecting the scalar potential 
of H, as in section 3.2, we perturb its magnetic part so that the static initial 
magnetic field Bo becomes the rotating magnetic field [^2) : 

B(t) = (BcosLut,Bsmut,B ). 

This system is closer to the systems used by other people to analyse the geo- 
metric phase, so it is important to determine its cyclic states. The Hamiltonian 
can be written 

In order to eliminate the time dependence, let us make the 'transition to the 
rotating frame' Q|, i.e., express the evolution operator as follows: 

u(t) = e -^ tL - e - iGt , 

where G is the time-independent Floquet generator: 

G = H(0) - loL z . 

As G is quadratic in v T = (r, p), the motion is determined by the kind of linear 
transformation induced on v by e~ %Gt : 

v(t) = e' iGt ve~ lGt = e At v, [iG, v] = Av. 
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This depends on the roots of the characteristic polynomial of A, Det(A7 — A), 
which become dependent of three dimensionless parameters 

\e\B \e\B uj c u 



«o 



2mcuj 2mcuj 2u> uj 

Our interest is centered in the case when all the roots of Dct(A7 — A) arc purely 
imaginary because in such a case G will induce 'confined' motions; this just 
restricts the parameters a, ao, w to some region in a — a — w space. By 
simplicity, we assume that the static fields of the Penning trap induce the PL 
with period r = 2T, i.e., 

2u) c 2oj c 4ao 

With this assumption, we can illustrate the classification of the 2-dimensional 
parameter space a — ao according to the nature of the motion induced on the 
charged particle (see figure 4) . The regions for which the motion is trapped are 
labelled as Tj, i = 1, . . . ,4, while the rest of the parameter domain produces 
deconfined motion 1231. 




Figure 4: classification of the a — ao plane into regions in which the charged particle 
performs confined motion (the regions labelled by Ti, i = 1,...,4) and deconfined 
motion (the rest of the plane) . 

From now on, let us suppose that the parameters a and ao belong to one of 
the regions Ti, i = 1, . . . ,4. In such a case, the Floquet generator G becomes 
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the superposition of three harmonic oscillator Hamiltonians: 

3 

i=i 

[A i ,Al}=e j S ij , a = ±l. 

Notice that some of the Ci,i = 1,2,3 could be negative, and in such a case 
there will be a global minus for the oscillator Hamiltonian contributing to G; 
this will be reflected in the spectrum of G which will not be bounded by below 
(see also §§]). 

Denote the eigenstates of G by \£ nin2n3 ), i-e. 

G\£ ni n 2 n 3 ) £ nin2 n 3 \£ nin2 n 3 ) . 

Let us assume that 

1-0(0)) = |£ ni „ 2 „ 3 ). 
The time evolution of this state is very simple: 

|^(t)) = e- iutL *e- iGt \S nin2n3 ) = e-^^e-^lSn^n,). 

Notice that this state is cyclic with period r = 2ir/ui 

|^(r)) = e~ l£ ^^ T \^{Q)). 

Its corresponding geometric phase is easily evaluated 

0ni n 2 n3 ^^{£nin 2 n 3 \L^> \£nin 2 n 3 } ■ 

It is possible to express the previous phase in terms of £„ in2n3 ; indeed, we have 
gotten a beautiful expression [^l], : 

R - _9 —f - _9 V" ( -l l \ du)i 

Pnin 2 n 3 ^ 7I " &ni n 2 n 3 — ^ / y I n i T ^ I ■ 

i— 1 ^ ' 

Thus, if we change slightly the rotation frequency u> of the field there will be 
a slight change in the levels of the Floquet generator G. The difference of the 
new level position and the old one at first order in u> is essentially the geometric 
phase. 



14 



David J. Fernandez C. 



6 Conclusions 

• We have seen that the dynamical manipulation, a procedure which looks 
for the variety of the operations that can be dynamically induced on 
a physical system, provides a more global notion of dynamics than the 
standard one. 

• We have seen also that the dynamical manipulation leads in a natural 
way to the study of some fundamental and practical problems in QM, 
as the squeezing transformation and the wavepacket reduction in a non- 
demolishing arrangement. 

• Finally, we have seen that those special dynamical processes called evo- 
lution loops provide the tools to easily evaluate the geometric phases for 
the corresponding cyclic evolutions or to simplify their calculation. 
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